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Intersubband response in a superlattice subjected to a homogeneous electric field (biased super- 
lattice with equipopulated levels) is studied within the tight-binding approximation, taking into 
account the interplay between homogeneous and inhomogeneous mechanisms of broadening. The 
complex dielectric permittivity is calculated beyond the Born approximation for a wide spectral 
region and a low-frequency enhancement of the response is found. A detectable gain below the 
resonance is obtained for the low-doped Ga^ls-based biased superlattice in the THz spectral region. 
Conditions of the stimulated emission regime for metallic and dielectric waveguide structures are 
discussed. The appearence of a localized THz mode due to BSL placed at the interface vacuum- 
dielectric is described. 

PACS numbers: 73.21. Cd, 78.45.+h, 78.67.-n 



a 

£ 

O 

(N 
> 
m 
O 
in 

o 

o 

-i— > 

B 

i 

-a 
c 

o 
o 



X 



I. INTRODUCTION 



The mechanisms of the stimulated emission due to in- 
tersubband transitions of electrons in different tunnel- 
coupled structures (monopolar laser effect) have been 
investigated during the previous decade (see Refs. in 
P, [H). As a result, both mid-IR and THz lasers via- 
bility has been demonstrated with the use of the scheme 
based on the vertical transport through quantum cascade 
structures, which incorporate injector and active regions. 
Since a population inversion appears in the active region, 
the stimulated emission occurs for the mode propagated 
along mid-IR or THz waveguide. 

In the case of a biased superlattice (BSL), the vertical 
current through the Wannier-Stark ladder, which takes 
place under the condition 2T <C e B (here e B /h is the 
Bloch frequency and T stands for the tunneling matrix 
element between adjacent quantum wells (QWs) 0,0]), 
does not change level populations. Thus, the consider- 
ation based on the golden rule approach gives zero ab- 
sorption. In contrast to this, for the wide miniband BSL, 
with the bandwidth 2T ^S> e B , a negative differential con- 
ductivity, i.e. gain due to Bloch oscillations, was stud- 
ied theoretically starting the 70s [5[ (see further results 
and references in Q) and demonstrated in recent exper- 
iments @, H- At the same time, a similar behavior of 
the THz response, including a crossover from gain to ab- 
sorption regime with detuning energy shifted through the 
resonance, was reported for the BSL with tight-binding 
electronic states [9(. This contradiction with the sim- 
ple quantum picture, where a zero response should take 
place, and the question about THz gain without inversion 
beyond the Born approximation are discussed in [lol . [ll[ . 
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FIG. 1: Transitions between broaded levels in BSL of period Z 
with the Wannie-Stark ladder (a). Peak and dispersive con- 
tributions of these transitions to the dielectric permittivity, 
Ae„ (solid and dashed curves are correspondent to i?e(Ae^) 
and Im(Ae^) respectively) (6). Geometries of the THz waveg- 
uides with the BSL placed between the ideal metallic mirrors 
(c), and at the interface vacuum-dielectric (d). 



To the best of our knowledge, the consideration of the 
response of a tight-binding BSL is not performed yet in 
spite of both real and imaginary contributions to the di- 
electric permittivity, Ae w , appear to be essential for wide 
miniband BSLs 

In this paper, we evaluate the response of a BSL placed 
on a high-frequency electric field taking into account both 
homogeneous and inhomogeneous mechanisms of broad- 
ening exactly. Within the tight-binding approach, which 
corresponds to the sequential tunneling picture (Fig. la), 
we analyze the frequency dispersion of complex dielec- 
tric permittivity (Fig. 16). The Green's function for- 
malism is used to describe both homogeneous and in- 
homogeneous mechanisms of broadening, and the quasi- 
equilibrium distribution of electrons over tunnel- coupled 
wells. We demonstrate the low-frequency enhancement of 
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the response under consideration. Further, the propaga- 
tion of the transverse magnetic (TM) mode along THz 
waveguides with the BSL placed between ideal metallic 
mirrors (Fig. lc) and at the interface vacuum-dielectric 
(Fig. Id) is considered. The appearance of a localized 
mode is founded for the second case. We also discuss the 
conditions for the stimulated emission regime for the two 
cases under consideration. 

The paper is organized as follows. In the next section 
we consider the THz response of BSL and in Sec. Ill we 
discuss the spectral dependencies of dielectric permittiv- 
ity, while in Sec. IV we consider the mode propagation 
for the BSL placed in the THz waveguide. The last sec- 
tion includes a discussion of the approximations used and 
conclusions. 



II. BASIC EQUATIONS 

Within the framework of the tight-binding approach 
we describe the electronic states in BSL using the in- 
plane inhomogeneous matrix Hamiltonian h rr i and the 
non-diagonal perturbation matrix due to a transverse 
field [E ± exp(— iwi) + c.c] written as [6h rr > exp(— iuii) + 
H.c.}: 



P 

2m 



Vr~. 



re. 



S rr > + T(6 rr '-i + <5 rr '+i), 



5h r 



ev ± 



E ± (r - r')((5 rr /_i + 5 rr > +1 ). 



(1) 



Here r — 0,±1, . . . stands as a QW number, p 2 /2m is 
the in-plane kinetic energy operator, and m is the ef- 
fective mass. The random potential energy of the r-th 
QW, V rx , is statistically independent in each QW and 
includes both short- and long-scale parts of potential. 
The Bloch energy, e B — \e\FZ, appears in (1) due to 
the shift of levels in the SL with period Z under a ho- 
mogeneous electric field F, and v ± = TZ/h stands for 
the transverse velocity [12j. The high-frequency current 
density, [/ w exp(— iu>t) + c.c], which involves oc lu^ 1 re- 
sponse and the intersubband contribution induced by the 
perturbation Sh rr i, is determined by: 



.e 2 n 2ev ± 

Iu = i E ± + l—r- 

muj L a 



^2 s J?\\( S Pr+lr ~ 5 Pr-lr 



where n is the electron concentration, factor 2 is due to 
spin, spy ... is the trace over in-plane motion, ((. . .)) is 
the averaging over random potentials V rx , and L is the 
normalization volume. 

The high-frequency contribution to the density matrix 
in Eq. (2), [Sp rr , exp(-iujt) + H.c], is governed by the 

independent linearized equations for 5p r = 5p r ± lr , see 
MM: 



. f (±> l ,i (±) f~ <±)f , 

-lLO0p r +^(/l r ±l<V> r ~ °Pr "r) 



±i-—F ± (p r± i 
nuj 



Pr) 



(3) 



with the in-plane Hamiltonian of the r-th QW written 
in the form h r = p 2 /2m + V rx + re B . Here we restrict 
ourselves to the consideration of only oc T 2 contribu- 
tions and the steady-state density matrix is written as 
p r . Next, we describe the electron states 
in the r-th QW by the use of the eigenstate problem 
(p 2 /2m + y rx )^ x = e ru ip^, where the quantum number 
v marks an in-plane state. Using this basis, we transform 
Eq. (3) into the form: 



(Sriln - £ru> ±E B -TlUJ - iXjSp 1 ^ (u, v') 



(4) 



Here A — > +0 and we apply the quasi-equilibrium distri- 
bution of r— th QW, p r — fp2/2m+v rx i where f e is the 
Fermi function with chemical potentials p, and temper- 
atures T e , which are identical for any QW. 
The current density (2) is given by 



.e 2 n 2ev ± 

4, ^ i E ± + i— - 

moj L 6 



(5) 



and the transverse conductivity, a^, introduced accord- 
ing to the standard formula I u = o~^E ± , takes the form: 



e 2 n 2(ev ± )' 2 

1 +l tT~ 

muj ujL a 




- fs r „,) (6) 
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where Q" u r * = J dx^ x *^, x is the overlap factor. We 

have replaced r — > r + 1 , v <-> v' in the second addendum. 
In order to check the non-singular behavior of o~ u at to — » 
0, one has to utilize the relation 



n 
m 



L 3 **r+l,r £ 



fe 



r+li/ 



= 0, (7) 



which can be proofed with the use of the definition 
n = 2V„ / Er „ and the relation v 2 ± Q^, = \{rv\v z \r' v')\ 2 , 
see [14]]. Using (7) we obtain the transverse dielectric 
permittivity, = e + iAira^/oj, in the form: 



2n(2ev ± ) 




3 



and intersubband transitions give a finite contribution to 
(8) in the static limit. Finally, we transform [l5| the de- 
nominators in (8) and rewrite Ae^ = — e through 
the spectral density function [ID, [lfij], .4 r>e (x, x') = 
J2v ' ( / ; rxV'rx*^( £ — £ rv), m the following form 



Aet = 



2 7 r(2ev_ L ) 5 

LU 2 L 3 



ds 



de'{f e -f e ,) 

-(e - e' - As- - iXy 1 - (e - e' + Ae+ + iA)" 1 ] 

x J dx J <bd ^((A+i, e (x,x')A^(x',x))) (9) 

with As± = Huj ± e B . Thus, we have evaluated the ex- 
pression for the transverse response taking into account 
the scattering processes exactly. 

Further, we perform the averaging over short-range 
and large-scale potentials in the last factor, keeping in 
mind that the averaged characteristics of scattering pro- 
cesses, both for homogeneous and inhomogeneous mecha- 
nisms, do not dependent on the QW number r. Following 
Eqs. (11 - 13) of Ref. [HI we obtain 



^£<<A + l,e(x,x')A,e<(x',x))) 



pOO 

g y o dSA(Z-e)A(Z-e f ), 



(10) 



where the averaged spectral density, A(E), is written in 
the integral form: 

Here T = y/Jw^j is the inhomogeneous broadening en- 
ergy due to the large-scale part of the potential in the 
r-th QW, uy x , and 7 stands for the homogeneous broad- 
ening energy. In (11) we consider the case of scattering 
by zero-radius centers when 7 does not depend on e, p or 
x and the shift of levels, which is logarithmically diver- 
gent without a small-distance cutoff [171 ] , is included into 
the zero point energy. The simple analytical expressions 
for the spectral density peak take the form in the limiting 
cases: 



A(E) 



7 /[7r(£ 2 + 7 2 )], 
exp(- J B 2 /2r 2 )/(\^ : r) 



r = 

7 = 



(12) 



and A(E) transforms from a Lorentzian towards a Gaus- 
sian line shape upon an increase in the contribution of 
the inhomogeneous broadening. 

According to Eqs. (9, 10), one needs to consider the 
convolution of spectral densities, D{e,e') = J °° d£A(£ — 
e)A(^ — £•'). For the collisionless case, when (12) is re- 
placed by the (5-function, one should replace D(e,e') by 
6{e + e')5(e — e'). Modifications of D{e, e') due to broad- 
ening are shown in Fig. 2 for the cases 7 = T and 7 = 3r. 
Using the variables e — e' and (s+s')/2 we obtain D(e, e') 
as a peak of width of the order of broadening with respect 
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FIG. 2: (Color online) Dimensionless convolution of spectral 
densities, 7-D(e, e'), plotted for 7 = V (a), and 7 = 3 F (b). 



to e — e' and D(e, e') is suppressed at negative (e + e') /2 
[HI . Note that tails of peaks are suppressed if the inho- 
mogeneous contribution increases. 



III. FREQUENCY DISPERSION 

Here we consider the frequency dispersion of dielectric 
permittivity tensor e|]; x . The in-plane component of per- 
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mittivity, e^, is written through the Drude conductivity 
°w = je 2 n/w(w + ij/h) in the form: 

e "= e r?^ ^ 

muj{u) + ij/n) 

with the homogeneous relaxation frequency, j/h, intro- 
duced in Sec. II. The transverse component is given by 
Eq. (9), which can be written through D(e,e') and the 
2D density of states, p 2D , as follows: 



± n(2ev ± ) 2 



de 



de' 



xD(s,s')(f £l -f e )[2/(e-s' + e B ) (14) 
-(e - e' - Ae- - iX)- 1 - (e - e' + Ae+ + iX)- 1 ] . 

The Fermi distribution, / e , is connected to the averaged 
concentration, n 2D = nZ ', according to the standard re- 
lation n 2D = f™^ dejepe = p 2D def e / °° d£A(£, - e), 
where p e is written through A{E) given by Eq. (11). 

We turn to estimates of the dielectric permittivity for 
the GaAs/Al , 3 Ga Q . 7 As BSL with a period Z = 170 A 
and with a tunneling matrix element T — 2 meV, which 
corresponds to a barrier width of 45 A. The level split- 
ting energy e B = 10 meV corresponds to a transverse 
electric field F = 5.9 kV/cm. Figure 3 shows the spectra 
of Ae w for the cases: 7 = T = 0.8 meV, and 7 = 3r with 
r = 0.4 meV. We have used in calculations a chemical 
potential p — 1.5 meV and temperatures T e — 0.5 meV 
and 4.5 meV. The corresponding 2D electron densities 
are n 2D = 6.1 x 10 10 cm" 2 (7 = T and T e = 0.5 meV), 
n 2D = 6.8 x 10 10 cm' 2 (7 = 3r and T e = 0.5 meV), and 
n 2D = 1.2 x 10 11 cm~ 2 (T e = 4.5 meV, both broaden- 
ing cases). Both the real (Fig. 3a) and imaginary (Fig. 
36) parts of Ae^ show a decrease of their peak value 
with increasing temperature. A non-symmetric shape 
of Re (Ae w ) appears due to the singular (cx oj~ 2 ) factor 
in Eq.(14). As a result, Im(Ae ul ) increases in the low- 
energy region. Homogeneous broadening, 7, has a bigger 
influence in the width of the permittivity dispersion than 
the inhomogeneous one, T. 



IV. ELECTRODYNAMICS 

Next, using the contribution to dielectric permittivity 
discussed above, we consider the in-plane propagation 
of a THz mode localized at the BSL, as it is shown in 
Figs. l(c, d). Since the only z-polarized component of 
the field is coupled to the BSL, we consider the TM- 
mode propagating along the OX-direction with the non- 
zero components [E n (z)e lku ' x : 0, E ± (z)e lkujX ]. The wave 
equation for these components can be transformed into 
the system (l9j : 





d 2 - 




' dz 2 _ 







E ± (z) 



= ik. 



— ik u 



dE ± (z) 

dz 
dE,{z) 

dz 



(15) 
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FIG. 3: (Color online) Frequency dispersion of the real (a) 
imaginary (b) parts of dielectric permittivity Ae^ . For the 
case 7 = T = 0.8 meV one uses p — 1.5 meV, T e — 0.5 meV 
(solid line) and p = 1.5 meV, T e = 4.5 meV (dot-dashed line). 
For the case 7 = 1.2 meV, T = 0.4 meV one uses p = 1.5 meV, 
T e = 0.5 meV (dashed line) and p = 1.5 meV, T e = 4.5 meV 
(dotted line). 



where e\£ is the dielectric permittivity tensor of the 
layered media with e\ u given by Eq. (13) and e^w = 
e z + Ae^ given by Eq. (14). 

Using the relation E ± (z) = ik u [dE ]] (z)/dz] 
/[ejr(uj/c) 2 — k 2 ] one may write the closed wave 
equation for (z) in the form: 



dz e^(u>/c) 2 — k 2 dz 



E t (z)=0. (16) 



We consider below three-layer structures with the BSL 
placed at \z\ < d/2. From Eq. (15) one obtains the 
wave equations with constant coefficients complemented 
by boundary conditions at z ± d/2: 



dE,{z) 



dz 



2=±d/2+0 
z = ±d/2-0 



(17) 
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and the continuity conditions: E l{ ( z )\ Z z= ±^/2-o = ^- 
In addition, the problem should be complemented by 
boundary conditions at \z\ ^> d. Below we restrict ourself 
by the cases of an ideal metallic waveguide and a THz 
mode localized at the interface vacuum-dielectric. 



A. Metallic Waveguide 

For the ideal metallic waveguide of width L we involve 
the additional boundary conditions E ]t (z = ±L/2) = 
and the wave equation (16) takes the form: 

(■^ + kAe 11 (z) = 0, \z\>i (18) 
(J^ + n 2 ^j E {l (z) = 0, \z\ < ^ 

where k and k± are determined from k 2 — e{uj /c) 2 — k 2 
and k 2 = e"(w/c) 2 — k 2 e^/e ± , respectively. Here e h± are 
the components of the BSL dielectric permittivity and 
e is the dielectric permittivity of the media inside the 
waveguide. 

We search the asymmetric solution of Eq. (17), which 
corresponds to the symmetric transverse field E ± (z), in 
the following form: 

)E W sin k (-| - z) , | < z < § 
£ sl sin Ki z, |z|<| , (19) 

-£' lu sinK(-| + z) , -|>2>-^ 

where the coefficients E w and E SL are determined from 
the above boundary conditions. The solvability condition 
gives us the dispersion relation: 

e L — d n ± d e" L — d n ± d 

— cos k — - — sin — - — I sin n — - — cos — — = 0, 

f% z z hb j_ z z 

which determines the complex longitudinal wave vector 
koj for the given BSL parameters and w. 

We solve the dispersion equation (20) using e 1,1 cal- 
culated in Sec. Ill for the BSL width d — 4/im and 
L = 40/im. Considering the lowest mode when the real 
part of k u appears to be close to y/eui / c, we have plotted 
Im(k LU ) in Fig. 4 for the above cases (7 = T and 7 = 3r). 
There are two perfectly defined spectral regions: a gain 
region with Im{k u ) < for huj < e B , and a damping one 
Im(k u! ) > for huj > e B . A marked increase of the gain 
can be seen when temperature T e decreases. 



B. Localized Mode 

For the case of a BSL placed near the vacuum-dielectric 
interface, we use Eqs. (15, 16) with the additional bound- 
ary conditions E^{z — > ±00) = 0. The wave equations 
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FIG. 4: (Color online) Frequency dispersion of Im(k^) for a 
metallic waveguide. Curves are marked the same as in Fig. 3. 

for vacuum, BSL and substrate regions take the forms: 

(£-«j)nw-o, .>£ 

(s? + "") i5 " w = °' w< f' (21) 

(£-«:)*<.) -0, .<-£ 

where k± is introduced in Sec. IIIA, k 2 — (oj/c) 2 — k 2 
and k 2 = e s (w/c) 2 — k 2 with the dielectric permittivity 
of substrate, e s . The field distribution is given by 

( E + e~ K ^ z - d / 2 \ z>d/2 
E^z) = I e+e lK ^ z + e-e- tK ^ z , \z\ < f , (22) 
[ E_eM z+d ' 2 \ z<~d/2 

where Re(n StV ) should be positive. The dispersion rela- 
tion is obtained from Eqs. (21, 22) in the form: 











(-" 


\- — ] cos n ± d — ( 


K ± 







The solution of this dispersion equation is performed 
for the above listed (7,T) and (fi, T e ) pairs. We use 
also the BSL width d = 4/im and e s = 3.7, which cor- 
responds to SiC"2 substrate. Figure 5 shows Im{k UJ ) 
for the cases under consideration. As in the metallic 
case, two regions can be selected: a gain region with 
Imiku) < for Tiuj < e B , and a damping one where 
Im{k w ) > for Tuo > e B . Contrary to the metallic case, 
gain for the dielectric waveguide exists even for low en- 
ergy values, at Tilo < 5 meV. Since Re(K v ^ s ) >> Im(K VtS ), 
the transverse size of mode is determined by i?e(K„ iS ) _1 
which are varied from 11.3 /im at low- frequency region to 
-Re(K„) -1 =9.45 /Ltm and i?e(K s ) _1 =5.1 /im at hui = 20 
meV. These results are weakly dependent on tempera- 
tures and the broadening cases under study. 
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FIG. 5: (Color online) Frequency dispersion of Im(k bJ ) for a 
dielectric waveguide. Curves are marked the same as in Fig. 3. 

V. SUMMARY AND CONCLUSIONS 

We have considered here the intersubband response of 
a biased superlattice on a THz radiation field beyond 
the Born approximation. Taking into account the inter- 
play between homogeneous and inhomogeneous broaden- 
ing we have analyzed the spectral and temperature de- 
pendencies of the complex dielectric permittivity in the 
low-doped BSL. We have found the low-frequency en- 
hancement of the dispersion of complex dielectric per- 
mittivity and have estimated the conditions to obtain 
the stimulated emission regime. The enhancement of the 
emission due to the THz waveguide effect is also consid- 
ered for the cases of the BSL placed between ideal metal- 
lic mirrors and at a vacuum-dielectric interface. The ap- 
pearence of the localized THz mode due to BSL placed 
at the interface vacuum-dielectric is described. 

Let us briefly discuss the assumptions used in the 
present calculations. The main restriction of the results 
is the description of the response in the framework of the 
tight-binding approach (within an accuracy of the order 
T 2 ) which is valid under the condition e B > 2T and is 
satisfied for the numerical estimates performed. Note 
that beyond the Born approximation the broadening can 



be comparable to the averaged electron energy. Next, in 
spite of the general equations (10) and (14) are written 
through an arbitrary spectral density function, with the 
use of statistically independent random potentials in each 
QW, final calculations were performed for a model that 
includes scattering by zero-radius centers and large-scale 
potential. Such a model describes the interplay between 
homogeneous and inhomogeneous mechanisms of broad- 
ening . Other approximations we have made are rather 
standard. We restrict ourselves to the case of uniform 
biased field and QW population neglecting a possible do- 
main formation caused by the negative differential con- 
ductivity of BSL 

One can avoid instabilities in a 
short enough BSL because the THz modes propagate in 
the in-plane directions. The Coulomb correlations, which 
modify the response as the concentration increases, are 
not taken into account here. This contribution, as well 
as the consideration of an intermediate-scale potential, 
requires a special attention in analogy with the case of 
a single QW [2l[. Finally, the simplified description of 
the ideal (without any damping) waveguide structure is 
enough to estimate the characteristic planar size of a de- 
vice suitable for THz stimulated emission: one have to 
compare the maximal negative value of 7m(/c w ) obtained 
with a damping length calculated for similar waveguides, 
see [22|]. Note, that more complicate waveguides (see 
Refs. in (23|) may be effective 

To conclude, the simplifications listed do not change 
the peculiarities of the THz response or the numerical 
estimates given in Sec. IV. It seems likely that the con- 
tribution of Re(Ae UJ ) can be found experimentally. More 
detailed numerical simulations are necessary in order to 
estimate a potential for applications of BSL as a THz 
emitter. 
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